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ABSTRACT 


The  dark  current,  impulse  and  square-pulse  response 
measurements  of  photoconductive  devices  fabricated  from  two 
different  types   of  materials.  Gallium  Arsenide  with  Chromium 
dopant  (GaAs:Cr)  and  Indium  Phosphide  with  Iron  dopant 
(InP:Fe)  are  reported.   These  devices  have  been  subjected  to 
irradiation  from  the  Naval  Postgraduate  School  S-band 
Electron  Linear  Accelerator  (LINAC)  with  an  energy  of  100  MeV 

at  room  temperature.   Fluence  ranged  between  10    and  10 

2 
electrons/cm  .   Dark  current  decreases  with  increasing  fluence 

for  the  GaAs:Cr  devices  whereas  InP:Fe  shows  an  increase  in 

the  dark  current.   Both  types  of  materials  exhibit  extremely 

fast  impulse  response  after  the  irradiation.   Electron 

mobility,  drift  velocity  and  response  speed  decrease  with 

increasing  fluence.   Response  speeds  <  100-ps  are  achieved 

by  fast  carrier  relaxation  in  the  semiconductor  due  to  the 

introduction  of  trapping  and  recombination  centers  resulting 

from  the  irradiation  damage.   The  GaAs:Cr,  unlike  the  InP:Fe, 

more  closely  follows  the  longer  square-pulse  exhibiting  no 

nonlinearity .   All  results  are  consistent  with  previously 

investigated  neutron  irradiated  devices. 
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.  I.   IJ^TRODUCTION 

There  has  been  recent  interest  in  the  development  of 
bulk-semiconductor  photoconductors  with  picosecond  carrier 
relaxation  time  to  be  used  as  radiation  detectors.   Data 
on  the  response  of  semiconductor  devices  to  radiation 
dam.age  is  required  for  design  and  development  of  their  use 
in  a  radiation  environment.   These  environments  include 
space  applications,  accelerators,  nuclear  reactors  and 
nuclear  explosions,  each  of  which  involve  specific 
irradiating  particles.   Compact  radiation  detectors  with  a 
high  degree  of  sensitivity  to  photons,  x-rays,  soft  x-rays, 
gamma  rays,  electrons,  protons  and  neutrons  are  being 
developed  by  the  Los  Alamos  National  Laboratory  Electronics 
Division.   The  current  response  of  these  devices  is  propor- 
tional to  the  incident-radiation  intensity  {Ref.  IJ . 

Particles  passing  through  devices  deposit  part  of  their 
energy  into  ionization  and  the  remainder  into  displacement. 
The  interaction  with  matter  depends  on  properties  of  both 
the  particle  (mass,  charge,  kinetic  energy)  and  the  target 
(mass,  charge,  density).   Bulk  device  characteristics  are 
degraded  by  displacement  damage  which  can  decrease  carrier 
mobility,  carrier  concentration  and  carrier  recombination 
lifetime. 

In  the  development  of  these  bulk-semiconductor  radiation 
detectors,  irradiating  particles  are  being  utilized  to 
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damage  the  device  thereby  creating  trapping  and  recombination 
centers.   The  desired  characteristic,  specifically,  a 
decrease  in  carrier  relaxation  time  resulting  in  response 
time  <  100-ps  is  obtained  while  maintaining  the  sensitivity 

of  the  device.   In  particular,  this  development  entails  the 

13 
use  of  neutrons  with  a  thermal  neutron  flux  of  9.0  x  10 

2 
(neutrons/cm  -sec)  for  periods  ranging  from  2000-4000 

seconds.   Neutron  transmutation  doping  in  the  bulk- 
semiconductor  device  is  impeded  by  the  filtration  of  low- 
energy  neutrons  [Ref .  IJ .   This  source  is  provided  by  the 
OMEGA  WEST  REACTOR  located  at  the  Los  Alamos  National 
Laboratory . 

In  this  thesis,  discussion  is  presented  on  the  radiation 
effects  of  high  energy  electrons  on  two  different  types  of 
materials.   Before  irradiation,  measurements  of  dark  current, 
impulse  response  (source  response  shorter  than  detector 
carrier  relaxation  time)  and  square-pulse  response  (source 
response  longer  than  detector  carrier  relaxation  time)  were 
made.   The  same  measurements  were  made  after  the  irradia- 
tion to  determine  changes  in  carrier  relaxation  time  giving 
an  indication  of  the  resulting  damage.   These  measurements 
were  carried  out  utilizing  the  processes  of  photoconductivity 

The  purpose  of  the  experiment  reported  here  is  to  deter- 
mine the  effect  of  100  MeV  electron  irradiation  on  these 
bulk-semiconductor  radiation  detectors  and  the  degree  to 
which  the  displacement  effects  can  be  compared  with  those 

12 


of  neutron  irradiated  devices.  It  is  desirable  to  obtain 
similar  results  with  the  charged  particle  source  (100  MeV 
electrons),  as  has  been  demonstrated  with  the  neutron  source. 
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II.   THEORY 

A.   IRRADIATION 

1.   Interaction  of  Electrons  with  Matter 

Electrons  incident  on  a  semiconductor  interact  with  the 
atoms  primarily  through  the  process  of  Coulomb  scattering.   The 
electron  path  may  be  deflected  through  large  angles  and  ranges 
may  vary  widely  even  for  electrons  of  the  same  energy.   The 
electron  energy  loss  per  collision  may  be  very  high. 

At  an  energy  of  100  MeV,  the  electrons  lose  energy 
partly  due  to  ionization  of  the  atoms  and  partly  due  to 
Bremsstrahlung  (radiation  loss) .   Stopping  power  is  defined  as 
the  amount  of  energy  lost  by  a  particle  per  unit  length  of  path 
through  the  stopping  material.   For  electrons,  the  total  stopping 
power  is  divided  into  two  components:   1)  collision  stopping 
power — the  average  energy  loss  per  unit  pathlength  due  to  inelas- 
tic Coulomb  collisions  with  bound  atomic  electrons  of  the  medium 
resulting  in  ionization  and  excitation,  and  2)  radiative  stopping 
power — the  average  energy  loss  per  unit  pathlength  due  to  the 
emission  of  bremsstrahlung  in  the  electric  field  of  the  atomic 
nucleus  and  of  the  atomic  electrons  [Ref .  2] .   As  the  electron 
traverses  its  path,  energy  going  into  ionization  and  excitation 
is  absorbed  by  the  target  relatively  close  to  the  path.   Energy 
lost  through  bremsstrahlung  travels  farther  from  the  path  before 
being  absorbed. 
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a.   Electron  Mass  Collision  Stopping  Power 

The  electron  mass  collision  stopping  power  is 
calculated  according  to  Bethe ' s  stopping  power  theory  using  the 
formulation  of  Rohrlich  and  Carlson  [Ref.  3] 

1  ?)F  2ttN  r  mc  Z  ^  _ 

^[||]  =   ^ {ln(^)^  +  ln(l+J)  +  F  (T)  -  6} 

^  collision        3  A 

(1) 

where : 

2      t2 

F  (T)   =   (1-3  )  [1  +-g-  -  [2T+l]ln  2] 

The  collision  stopping  power  of  a  compound  can  be  approximated 
by  the  weighted  sum  of  stopping  powers  of  the  atomic  constitu- 
ents of  the  compound.   In  formula  1,  the  following  replacements 
are  applicable. 


Z/A  by  <Z/A>   =  \   w. (Z./A. ) 

J     -     ^ 


In  I   by   In  I   =   [  J]  w .  (Z  ./A  • )  In  I.]/<Z/A>  (2) 


6   by   6   =   [y  w. (Z ./A. ) 6 .]/<Z/A> 

^   J   :   3   D 


where: 


1)  w .  is  the  fraction  by  weight  of  the  j ' th  atomic 

constituent 

2)  Z.,  A./  I.  and  5.  pertain  to  the  j ' th  atomic 

constituent 
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Properties  of  the  medium  are: 
Z   =   atomic  number 

A   =   atomic  weight 

3 
p   =   density  (gm/cm  ) 

I   =  mean  excitation  energy  (eV) 

6   =   density  effect  correction 

Properties  of  the  electron  are: 

2 

T   =   T/mc   -  kinetic  energy  in  units  of  rest  mass 

3   =   v/c   -   1 
Other  parameters  are: 


2 

2ttN   r   mc        =       .  153536  ( t ) 

a  e  mole 


Listed  under  the  properties  of  the  medium  are  the 
mean  excitation  energy  and  the  density  effect  correction.   The 
mean  excitation  energy  is  the  logarithmic  average  of  the  exci- 
tation energies  of  the  medium  weighted  by  the  corresponding 
oscillation  strengths.   One  assumption  used  in  determining  the 
mean  excitation  energy  is  the  use  of  the  I-values  of  the  ele- 
mental substances  for  the  I-values  of  the  atomic  constituents 
of  the  compound.   This  neglects  the  molecular  binding  energies. 
The  density  effect  correction  takes  into  account  the  polarization 
of  the  medium  v/hich  screens  the  electric  field  acting  on  the 
relativistic  particle.   The  effect  is  a  reduction  of  collision 
stopping  power.   [Ref.  2] 
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b.   Electron  Mass  Radiative  Stopping  Pov/er 

The  electron  mass  radiative  stopping  power  can  be 
expressed  using  the  following  formula  [Ref.  2]. 


IrSE,  N  ar  EZ    „  t  (})^(t)   , 

^[t^J  =    a   e     .n       ^     1  ^     rad,       ,-.. 

P  ^^  radiative     A ^  ^^^rad^^  ^  Z  -jT-^ ^       ^^^ 

^  ^"^rad 


Properties  of  the  medium  (compound)  are: 

A   =   w.Al  +  W.A2 
Z   =   w. Zl  +  w. Z2 

where: 

1)  w.  is  the  fraction  by  weight  of  the  j ' th  atomic 

-^  constituent 

2)  Al   =   atomic  weight  of  element  1 

3)  A2   =   atomic  weight  of  element  2 

Properties  of  the  electron  are: 


2 

E   =   T  +  mc    =   100.511  MeV 


Other  parameters  are: 

e  n 

^    (t)   ,  and  $  (t)   ,  are  dim.ensionless ,  scaled,  radiative 

energy  loss  cross  sections 
N^   =   6.022045  X lo^^/mole 

r^   =   7.940775  x 10^^  cm^ 
e 

a    =   Z/137 
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Implicit  in  the  formulation  of  the  electron  mass 
radiative  stopping  power  is  electron-nucleus  bremsstrahlung 
(emission  of  a  photon  due  to  the  interaction  of  the  electron 
with  the  screened  Coulomb  field  of  the  atomic  nucleus)  and 
the  electron-electron  bremsstrahlung  (Coulomb  interaction  with 

one  of  the  atomic  electrons) .   The  ratio  of  the  scaled,  radia- 

e        n 
tive  energy-loss  cross  sections  (6  (t)   V*  (t)   J  is  assumed 
^^  rad       rad 

to  be  one  and  values  for  {^    (t)   -.)    are  approximated  from  least- 
squares  curves  fitted  from  theoretical  points  based  on  Davies, 
Bethe  and  Maximon  (1954)  and  Olsen  (1955)  .       [Ref .  2] 
c.   Energy  Deposition--Dose  Measurements 

The  range  of  the  electron  is  defined  as  the  average 
path  length  the  electron  travels  until  it  is  stopped  by  the 
medium  as  a  result  of  energy  loss.   The  electrons  are  assumed 
to  lose  energy  continuously.   If  the  range  is  divided  by  the 
density  of  the  material,  the  distance  an  electron  will  travel 
can  be  determined. 

Dose  represents  the  total  amount  of  energy  deposited 
to  the  medium.   It  is  given  in  rads  (100  ergs/gram)  and  must  be 
specified  for  the  particular  material.   The  front  surface  dose 
(thin  sample)  is  expressed  by  the  following  formula  [Ref.  4]. 


R   =   1.6  X 10    $[-[3^3  i  (4) 

^  collision 


where : 


2 

electron  fluence  (electrons/cm  )  is  determined 

from  formula  (24) 
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2.   Lattice  Defects 

The  type  of  lattice  structure  inherent  in  the  study  of 
these  devices  is  the  zinc-blende  structure.   In  the  diamond 
lattice  structure,  each  atom  lies  in  the  center  of  a  tetra- 
hedron formed  by  the  four  nearest  neighbors.   This  structure 
is  restricted  to  elements.   The  zinc-blende  structure  is  similar 
to  the  diamond  lattice  with  the  exception  that  the  two  nearest 
neighbors  are  occupied  by  different  elements. 

There  are  various  types  of  lattice  defects  that  can 
occur  from  100  MeV  electron  irradiation.   Considering  only  point 
defects,  the  lattice  defects  are  broken  up  into  two  main  cate- 
gories which  include  isolated  defects  and  cluster  defects. 
Within  isolated  defects,  the  simplest  type  of  defect  assumed  is 
the  Frenkel  defect.   This  is  an  empty  lattice  site  (vacancy) 
and  an  atom  occupying  an  interstial  position.   The  Frenkel  defect 
requires  an  amount  of  energy  given  to  the  atom  needed  to  displace 
it  to  the  metastable  intersitial  site.   The  displacement  process 
depends  on  the  energy  of  the  recoil  and  its  direction  with  respect 
to  other  lattice  atoms  [Ref .  5J .   Other  isolated  defects  include 
divacancies  which  result  when  a  high  energy  recoil  atom  imparts 
sufficient  energy  to  displace  a  secondary  atom  along  with  the 
primary  atom.   The  more  complicated  cluster  defects  involve 
regions  containing  large  numbers  of  isolated  defects. 

For  binary  compounds  (zinc-blende  structure) ,  the  number 
of  possible  lattice  defects  is  considerably  higher  than  for 
elements.   The  displacement  or  replacement  of  atoms  may  result 
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in  eight  types  of  point  defects:   two  types  of  vacancies,  four 
types  of  interstitial  atoms  (surrounded  by  similar  or  dissimilar 
atoms)  and  two  types  of  substitution  defects  [Ref.  6]. 

In  the  following  sections,  the  irradiation  induced 
processes  of  atomic  displacements  will  be  considered  on  a  micro- 
scopic level.   Both  single  atom  and  multiple  atom  displacements 
are  discussed  along  with  the  assumptions  made  in  determining 
the  total  number  of  displaced  atoms.   A  complete  review  of 
damage  theory  is  found  in  Seitz  and  Koehler  [Ref.  7] . 

a.   Atomic  Displacement  Cross  Sections 

For  100  MeV  incident  electrons  (relativistic) , 
enough  energy  is  transferred  to  the  lattice  atoms  to  knock  them 
from  their  normal  lattice  position  into  an  interstitial  site. 
These  primary  knock-on  atoms  may  in  turn  depart  some  of  their 
energy  to  form  secondary  displaced  atoms.   Radiation  damage 
theories  concentrate  on  these  few  elastic  collisions. 

In  the  relativistic  range,  McKinley  and  Feshbach 
[Ref.  8]  formulated  an  expression  for  lower  Z  elements  using 
the  Rutherford  differential  cross  section. 


da  ^   =   R  ^  da   ^,  (5) 

mf      mr    ruth 


where 

R^f   =   (1-3^)  [1-3^  ;^  +  Tra3{(^)^^^  -^}] 

m        mm 

The  differential  cross  section  of  the  atom  must  be  integrated 
over  all  transferred  energies.   Additionally,  the  assumption 
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is  made  that  recoil  energies  below  the  threshold  displacement 
energy  (T.)  have  a  probability  of  zero  for  displacing  an  atom 
whereas  those  recoil  energies  above  the  threshold  displacement 
energy  (T. )  have  a  probability  of  one.   For  thin  samples,  the 
energy  loss  is  approximately  constant  over  the  sample  thick- 
ness.  The  displacement  cross  section  is  dominated  by  small 
angle  scattering  at  100  MeV  energies  and  is  given  by 

2    T  T 

^mf   "   2.495x10    cm  -^-^LC— -l)-3  In  ^ 

3  Y    ^  d 

T  T 

+  TTa3{2[(;^)-^/^-l]  -In  ^}]  (6) 

d  d 

Parameters  are: 

2 

.,             2(E+2mc  )E   ,  „, 
T   =  maximum  recoil  energy   =  — ^ ^ — - —   (eV) 

E    =   energy  of  the  incident  electron  =  100  MeV 

m    =   rest  mass  of  the  electron 

M   =  mass  of  the  target  atom 

T    =   threshold  displacement  energy  (values  given  in 
°-  Table  1) 

a    =   Z/137 

If  the  primary  atom  has  sufficient  energy,  additional 
displaced  atoms  are  produced  forming  a  displacement  cascade. 
In  order  to  determine  the  effective  cross  section  for  both  pri- 
mary and  secondary  displacements,  v(T),  the  average  number  of 
displacements  per  primary  recoil  of  energy  T,  must  be  determined. 
A  number  of  assumptions  are  made.   [Refs.  5,9,10] 
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TABLE  1 
THRESHOLD  DISPLACEMENT  ENERGIES  [REF.  9] 

Substance GaAs InP    

Displaced  Atom     Ga        As         In         P 

Threshold 

Electron  Energy 

E,   (keV)       .     228        273        270         110 

D 

Threshold 

Displacement 

Energy  T^  (eV)      8.8        10.1       6.6         8.8 

Self-Dif fusion 

Energy  (eV)         5.6        10.2       3.85        5.65 


1.  A  step  junction  threshold  is  assumed. 

2.  Atom-atom  collisions  are  treated  as  hard-sphere  elastic 
interactions. 

3.  The  ordered  arrangement  of  atoms  in  the  lattice  structure 
is  not  considered. 

4.  Damage  is  considered  homogeneous. 

5.  Annealing  is  not  considered. 

6.  Glancing  collisions  are  not  considered  even  though  energy 
loss  may  be  incurred. 

7.  The  number  of  replacements  per  primary  are  not  considered. 
The  moving  atom  will  replace  the  struck  atom  if  the  latter 
receives  energies  greater  than  T,  and  the  former  retains 
energy  less  than  T . . 

8.  Lattice  atom^s  are  considered  at  rest. 

9.  Long  range  effects  from  other  atoms  are  not  considered. 
10.   There  is  no  account  for  ionization  losses. 

Kinchin  and  Pease  [Ref.  11]  consider  the  threshold 
energy  by  postulating  that  both  the  striking  and  struck  atoms 
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require  an  energy  T,  or  greater  after  the  collision.   This 
accounts  for  replacements  but  does  not  account  for  the  energy 
loss  of  the  primary  atom  in  climbing  out  of  its  potential 
well.   Figure  1  displays  the  assumptions  made  by  Kinchin  and 
Pease  in  calculating  the  average  number  of  displacements.   The 
Snyder  and  Neufeld  model  [Ref.  12]  does  not  account  for  replace- 
ment but  does  account  for  the  potential  binding  energy.   In 
this  model,  the  secondary  atom  loses  energy  T,  before  making 
subsequent  displacements.   The  average  number  of  displacements 
per  primary  recoil  of  energy  T.  is  given  by 


v(T)   =   1     T .  <  T  <  2T, 

d  —   —   d 


■prp 

=   7=-    2T   <  T  <  T. 
^d      d 


(7) 


Parameters  are: 

T    =  average  energy  transferred   =   T.  (In  T  /T,  -1  +  rra) 

T.   =  Ionization  Energy   =   ME-,/8m 

B    =  .5  (Kinchin  and  Pease);  .56  (Snyder  and  Neufeld) 

For  compounds,  the  ratio  of  the  atomic  masses 
influences  the  transferred  energy.   The  maximum  energy  trans- 
ferred is  A  times  the  energy  of  the  striking  atom  where  [Ref.  9] 

A   =   ^^^ 2  (S) 

(M^  +  M  ) 
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Figure  1 


Kinchin  and  Pease  assumptions  for  calculating 
the  average  number  of  displacements  [Ref.  10] 
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Figure  2.   Kinchin  and  Pease  assumptions  for  calculating 
the  average  number  of  replacements  [Ref.  10] 
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The  Kinchin  and  Pease  model  holds  for  compounds  pro- 
vided there  is  equal  probability  of  collision  for  the  secondary 
atoms.   Replacement  collisions  involving  unlike  atoms  are 
called  misplacements  and  account  for  a  substantial  portion  of 
the  damage. 

The  effective  cross  section  for  primary  and  secon- 
dary displaced  atoms  in  compounds  is  determined  by  multiplying 
(6)  ,     (7)  and  (8)  , 


^d   =   %f  ^(T)  ^  (^) 


b.   Number  of  Displaced  Atoms 

In  the  previous  section,  the  total  cross  section 
for  both  primary  and  secondary  displacements  was  determined. 

From  this,  the  density  of  Frenkel  defects  (for  thin  samples) 

3 
produced  per  cm   can  be  determined  given  the  incident  electron 


energy  and  fluence.   [Ref.  13] 


Np   =   *e  °d  "o  <"' 


Parameters  are : 

3 

Np   =   number  of  Frenkel  defects  per  cm 

Oj   =   total  displacement  cross  section 


$   =   electron  fluence 
e 


-.  p  N  #  atoms 

-3        a 


N_   =   number  of  lattice  atoms  per  cm 
o  ^  A 
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where 


3 

p    =   density  (gm/cm  ) 

23 
N^   =   Avogadro ' s  number  =  6.0  2  x lo    (atoms/mole) 
a 

A   =   atomic  number  of  the  target  molecule 


The  number  of  Frenkel  defects  determined  experi- 
mentally from  various  sources  is  always  less  than  the  number 
predicted  by  theory.   This  is  a  function  of  the  assumptions  made 
in  formulating  the  various  models, 
c.   Replacement  Collisions 

Replacement  collisions  occur  when  there  is  an  ex- 
change of  roles  between  the  primary  and  secondary  knock-on 
atoms.   For  binary  compounds,  the  exchange  can  be  between 
similar  or  dissimilar  atoms  and  will  occur  when  the  moving 
atom  imparts  sufficient  energy  to  the  stationary  atom  to  cause 
displacement  but  does  not  retain  sufficient  energy  to  escape 
the  potential  well  of  the  vacancy,  thereby  replacing  the  dis- 
placed atom. 

Dienes  and  Vineyard  [Ref.  10]  have  determined  the 
number  of  replacements  per  primary  atom  of  energy  T  to  be 


T 
y(T)   =   2^[1.614  In(^)  +  IJ    T  >  2T^  (11) 

d  r 


The  replacement  threshold  T   has  not  been  determined  experimen- 
tally, however.  Kinchin  and  Pease  have  suggested  a  ratio  of 
T./T   =10.   A  summary  of  the  Kinchin  and  Pease  assumptions  for 
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determining  the  average  number  of  replacements  is  found  in 
Figure  2. 

d.   Radiation  Effects 

The  macroscopic  electrical  properties  of  semicon- 
ductors are  affected  by  the  microscopic  defects  resulting  from 
electron  irradiation.   The  specific  type  of  damage  is  a  function 
of  many  different  variables  including  temperature  during 
irradiation,  energy  spectra,  impurity  and  dopant  concentrations, 
rate  and  history  of  the  irradiation  and  minority  carrier  injec- 
tion ratio.   The  introduction  of  energy  levels  within  the 
forbidden  gap  due  to  point  defects  will  be  considered. 

Point  defects  primarily  affect  the  electrical 
properties  by  their  ability  to  either  produce  or  remove  charge 
carriers,  change  carrier  mobility  or  both.   Several  models  have 
been  proposed  to  explain  the  behavior  of  Frenkel  pairs  in  semi- 
conductors, especially  in  silicon  and  germanium.   James  and  Lark- 
Horovitz  [Ref.  14]  have  suggested  that  interstials  act  as  donors 
below  the  conduction  band  with  two  donor  levels  assigned  whereas 
empty  acceptor  states  above  the  valence  band  are  produced  by 
vacancies.   The  relative  location  in  the  forbidden  gap  will 
ultimately  determine  the  net  effect  of  these  localized  states 
with  electrons  redistributing  themselves  to  fill  the  low  energy 
states.   Blount  [Ref.  15]  based  his  arguments  on  atomic  and 
molecular  orbitals  stating  that  each  defect  can  act  as  either 
donor  or  acceptor,  depending  on  the  position  of  the  Fermi  level. 
As  the  concentration  of  irradiation  defects  increase  to  a  level 
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which  is  considerably  greater  than  the  concentration  of  impuri- 
ties, the  Fermi  level  approaches  saturation  and  is  determined 
by  the  distribution  of  defect  levels  in  the  forbidden  gap. 
Table  2  lists  energy  levels  that  have  been  experimentally 
determined  for  InP  and  GaAs  [Refs.  9,16]. 

TABLE  2 
InP  AND  GaAs  ENERGY  LEVELS 

Band  Gap  Forbidden  Gap  Energy  Levels  in  eV  at  300 °K 

E„/  opt  (300°K)     with  Electron  or  Gamma  Irradiation 

G 

Below  Conduction  Band     Above  Valence  Band 
InP       GaAs  InP         GaAs         InP       GaAs 

1.260     lc430  0.285       0.13  0.35 


InP 

GaAs 

0.285 

0.13 

0.38 

0.52 

The  effect  on  the  processes  of  photoconductivity 
due  to  the  introduction  of  energy  levels  in  the  forbidden  gap 
during  electron  irradiation  will  be  discussed  in  the  following 
sections . 

B .   PHOTOCONDUCTORS 

1 .   Photoconductivity 

A  photoconductor  consists  of  a  bulk-semiconductor 
material  with  ohmic  contacts  (Figure  3) .   The  electrical 
conductivity  is  equal  to  the  product  of  the  number  of  charge 
carriers,  their  charge  and  their  mobility.   The  steady-state 
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Figure  3.   Picosecond  Photoconductive  Radiation  Detector 
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CONTACT 

\            TEMPORARY               / 
\C0NDUCTIVE  REGION   / 

SEMICONDUCTOR  MATERIAL 

Method  of  Operation.   The  region  between  the  ohmic 
contacts  is  illuminated  with  the  photon  source. 
The  photons  are  absorbed  creating  free  charged 
carriers  (electrons  and  holes)  which  produce  a 
temporary  conductive  path  between  the  contacts. 
The  full-width  at  half  maximum  (FWHM)  impulse 
response  measurement  is  a  measure  of  the  photo- 
conductor  response  speed.   It  is  a  function  of 
the  carrier  recombination  lifetime  in  the  active 
region.   [Ref.  22] 
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contribution  to  the  electrical  conductivity  resulting  from 
photon  irradiation  is  defined  as  the  photoconductivity.   In 
nonhomogenous  semiconductors,  the  ideal  of  an  "effective  mobility" 
is  used  to  describe  low  conductivity  barriers  which  limit  the 
flow  of  current  through  the  material  [Ref.  6].   Steady-state 
photoconductivity  is  expressed  as 


Aa   =   q(nAu*  +  pAy*)  (12) 


When  excited  by  an  external  photon  source,  charged 
carriers  are  generated  by  either  band-to-band  transitions  (in- 
trinsic) or  by  transitions  involving  forbidden  gap  energy  levels 
(extrinsic) .   It  is  the  extrinsic  transitions  which  will  be  of 
interest  in  this  experiment.   Photoconductivity  takes  place  by 
the  absorption  of  photons  of  energy  equal  to  or  greater  than  the 
energy  separation  between  the  forbidden  gap  levels  and  the 
conductance  or  valence  band  levels.   The  conductance  of  these 
devices  is  proportional  to  the  intensity  of  the  incident  pho- 
ton source  and  takes  place  via  a  temporary  conductive  region 
between  the  ohmic  contacts  (Figure  4)  resulting  from  the  crea- 
tion of  free  electrons  and  holes.   This  nearly  intrinsic  region 
consists  of  a  high  concentration  of  forbidden  gap  energy  levels 
acting  as  recombination  centers.   For  transient  photoconduc- 
tivity, the  picosecond  response  time  depends  on  the  carrier 
relaxation  time  and  not  on  the  spacing  between  contacts.   This 
allows  for  a  large  bias  across  the  device  thereby  increasing 
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the  sensitivity.   Quasi-neutrality  and  ideal  ohmic  contacts  are 
assumed.   [Ref.  17] 

Desirable  material  properties  for  radiation  detec- 
tor applications  include  short  carrier  recombination  lifetime, 
high  carrier  drift  velocity,  high  resistivity  (low  dark  current) 
and  high  carrier  mobility  [Ref.  18].   The  efficiency  of  the 
photoconductor  is  measured  in  terms  of  its  photoconductive  gain, 
response  time  and  sensitivity.   It  has  been  concluded  that  the 
transient  recombination  is  a  bulk  controlled  process  for  InP:Fe 
[Ref.  18] . 

2.   Recombination  Centers  and  Carrier  Relaxation  Time 

The  time  that  a  charged  carrier  is  available  to  contribute 
to  the  conductivity  is  defined  as  the  free  lifetime.   For  excited 
electrons,  it  is  the  time  spent  in  the  conduction  band.   For 
excited  holes,  it  is  the  time  spent  in  the  valence  band.   Both 
minority  and  majority  carriers  contribute.   The  free  lifetime 
is  limited  by  three  different  processes:   1)  extraction  of  the 
carriers  from  the  photoconductor  without  replacement  as  a  result 
of  a  voltage  bias,  2)  absorption  into  a  trapping  center,  3)  recom- 
bination with  a  carrier  of  opposite  sign  at  a  recombination 
center.   The  distinction  between  a  trapping  center  and  a  recom- 
bination center  is  statistical  in  nature  and  is  delineated 
as  follows:   in  a  trapping  center,  a  captured  carrier  has  a 
greater  probability  of  thermal  re-excition  to  a  free  state  whereas 
a  recombination  center  provides  the  captured  carrier  with  a 
greater  probability  of  recombining  with  a  carrier  of  the  opposite 
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sign.   A  demarcation  line  represents  that  level  which  gives 
equal  probability  for  either  event  to  occur.   Figure  5  displays 
the  relative  position  between  the  Fermi  level  and  the  demarca- 
tion level.   Area  I  represents  those  levels  (electron  traps) 
in  which  the  occupation  depends  on  the  thermal  equilibrium 
between  the  conduction  band  and  the  level.   The  occupation  of 
the  levels  in  Area  IV  (hole  traps)  depends  on  the  thermal 
equilibrium  between  the  valence  band  and  the  level.   Due  to 
the  high  density  of  occupied  levels  in  Area  II,  recombination 
is  taking  place  even  though  the  conduction  band  is  still  in 
thermal  equilibrium  with  these  levels.   [Ref .  19] 

Recombination  takes  place  via  bound  states  in  the  forbidde 
gap.   These  bound  states  are  comprised  of  impurities,  crystal 
defects,  vacant  lattice  sites  and  interstial  atoms  [Ref.  20]. 
Before  irradiation,  the  energy  levels  in  the  forbidden  gap  are 
due  to  the  chemical  doping  of  the  sample  and  determine  the 
number  of  charge  carriers.   Irradiation  produces  defects  result- 
ing from  atomic  displacements  and  introduces  energy  levels 
into  the  forbidden  gap  of  semiconductors.   These  energy  levels 
are  often  deeper  in  the  gap  than  those  dominating  conductivity 
and  serve  as  recombination  centers  for  electrons  and  holes.   The 
process  of  recombination  at  the  local  centers  is  influenced  by 
the  carrier  population  of  the  centers.   The  capture  of  the 
minority  carrier  is  generally  the  limiting  step  in  the  recom- 
bination process.   There  are  four  processes  through  which  free 
carriers  interact  at  these  recombination  centers:   1)  electron 
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Figure  5.   Semiconductor  Fermi  levels  and  demarcation 
levels  [Ref.  19] 
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capture — the  process  in  which  an  electron  falls  from  the 
conduction  band  into  an  neutral  localized  state,  2)  electron 
emission — the  emission  of  an  electron  from  the  negatively 
charged  localized  state  into  the  conduction  band,  3)  hole 
capture--the  falling  of  an  electron  from  the  negatively  charged 
localized  state  into  the  valence  band,  4)  hole.  emission--the 
excitation  of  an  electron  from  the  valence  band  into  the  neutral 
localized  state.   These  processes  are  illustrated  in  Figure  6. 
The  net  rate  of  recombination  is  defined  by  the  following 
equation : 


U   = 
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where  a   and  a   are  the  electron  and  hole  capture  cross  sections, 
n       p  ^ 

v^,  is  the  carrier  thermal  velocity,  N_^  is  the  recombination 
th  -^  '   t 

center  density,  E   is  the  recombination  energy  level  and  E.  is 
the  intrinsic  Fermi  level. 

Figure  7  displays  the  relative  position  of  the  Fermi 
level  (E^) ,  the  recombination  center  (E  ) ,  center  of  the  band 
gap  (E.)  and  the  demarcation  level  (En).   The  minority  lifetime 
of  a  free  carrier  varies  depending  on  the  relative  positions 
of  these  levels.   All  six  cases  assume  the  density  of  the 
recombination  centers  to  be  greater  than  the  density  of  the 
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Figure  6.   Recombination  processes  [Ref.  19] 
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Figure  7 


Relative  positions  between  Fermi  level,  recom- 
bination center  level,  center  of  band  gap  and 
demarcation  level  in  a  semiconductor  [Ref.  19] 
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free  carriers.   Additionally,  the  recombination  processes  apply 
to  all  levels  lying  between  the  Fermi  level  and  the  demarcation 
level,  with  an  exponential  decrease  in  the  contribution  to 
centers  outside  the  bounds  of  these  levels.   The  Fermi  level 
and  demarcation  level  move  in  opposite  directions  but  at  an 
equal  rate. 

The  carrier  relaxation  time  is  due  to  the  trapping  and 
recombination  of  the  free  carriers  in  the  conductive  region. 
A  mechanism  for  electron  and  hole  capture  in  the  InP:Fe  samples 
suggests  that  an  acceptor  state  .61  eV  below  the  conduction 
band  results  from  the  Fe  dopant.   This  acceptor  state  can  exist 
in  either  of  two  charged  configurations.   Electrons  are  captured 
via  the  neutral  state.   The  negatively  charged  state  acts  as 
a  recombination  center  with  the  capture  of  holes.   The  relaxation 
time  is  characterized  by  the  following  equations.   [Ref.  1] 


8n(t)/3t   =   f(t)  -  r^^  (14) 


8p(t)/8t   =   f(t)  -  r^^  (15) 


9n^_(t)/at   =   r^^  -  r^^  (16) 


^CT   =   ^e^e^TO^^^^(t)  ^^^^ 


TV 


V^O^n^_it)p{t)  (18) 
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where  n(t)  and  p(t)  represent  the  electron  and  hole  concentra- 
tions, n^_(t)  and  n   (t)  are  the  concentration  of  charged  and 

o 
neutral  acceptor  states,  f(t)  is  the  generation  rate,  r^^p  is 

the  capture  rate  of  electrons  by  the  neutral  acceptor  states, 
r    is  the  capture  rate  of  free  holes  by  the  negatively  charged 
acceptor  state,  v  and  v,  are  the  thermal  velocities  of  elec- 
trons and  holes,  and  a   and  o,     are  the  capture  cross  section 
for  electrons  and  holes  at  the  neutral  and  negatively  charged 
acceptor  state. 

Transient  photoconductivity  involves  the  assignment  of 
two  characteristic  times  to  the  forbidden  gap  energy  levels  to 
describe  the  recombination  process.   The  transient  time  con- 
stants for  small  excess  carrier  densities  in  a  nondegenerate 
semiconductor  near  thermal  equilibrium  are  described  by  the 
following  formulas.   [Ref.  17] 


An(t)   =   A^  exp  i-t/T^  +   B^  exp  i-t/x^^)  (19) 


and 


Ap(t)   =   A  exp  (-t/Tj^  +  B  exp  (-t/Tj^j^)  ,  (20) 


where  the  initial  conditions  determine  the  coefficients  A  ,  B  , 

n'   n 

A   and  B  .   The  equal  capture  rate  for  electrons  and  holes  by 
a   recombination  center  requires  a  readjustment  of  charges. 
The  time  needed  for  this  readjustment  is  denoted  by  t    and 
is  identical  for  both  electrons  and  holes.   For  the  condition 
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when  N^„  <<  n.,  t»,„  (the  main  recombination  lifetime)  does 
TT      1    MR 

not  equal  the  steady-state  lifetime.   [Ref.  17] 
3 .   Impulse  Response 

Impulse  response  occurs  when  the  source  response  is 
shorter  than  the  detector  carrier  relaxation  time.   The  conduc- 
tance of  the' device  is  the  result  of  photon  absorption  and  in- 
creases linearly  with  the  energy  of  the  radiation  pulse.   This 
absorption  produces  free  electrons  and  holes  that  provides  a 
conductive  path  between  the  metallic  contacts  (as  illustrated 
in  Fig.  4).   The  response  time  is  a  property  of  the  material 
and  depends  inversely  on  the  uncharged  recombination  centers. 
The  100  MeV  electron  irradiation  damages  the  device  in  such  a 
way  as  to  create  trapping  and  recombination  centers  which  pro- 
duce faster  response  times  and  increased  sensitivity. 

Desirable  properties  for  detector  application  includes 
a  peak  current  response  which  is  proportional  to  the  light 
intensity  and  a  response  speed  which  is  independent  of  light 
intensity  [Ref.  21] .   The  peak  current  response  is  indicated 
by  the  following  formula.   [Ref.  22] 

(1-R  )  q  V   E 

I   =   ^ —  (21) 

S  w 

Parameters  are: 

I    =   observed  peak  current 

-19 
q    =   electronic  charge  =  (1.6  xio     C) 

V-,   =   electron  drift  velocity  =  y  V    ,/S 
d  J'    I-   appl 

E    =   Energy  deposited  into  detector  by  incident 
radiation  =  (580  ±20%  pJ) 
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S    =   contact  spacing   =  (.05  cm) 

W   =   Energy  to  create  an  electron-hole  pair  =  (5  eV) 

R   =   Reflectivity  =  (30  percent) 

The  assumptions  implicit  in  this  formula  are  an  incident 
radiation  pulse  which  is  much  shorter  than  the  carrier  relaxa- 
tion time  of  the  device,  the  impedance  of  the  device  is  always 
much  greater  than  that  of  the  50  ohm  load  and  current  resulting 
from  the  free  holes  is  ignored  because  of  the  much  lower  drift 
velocity. 

The  electron  mobility  in  the  radiation  detector  may  be 
approximated  by  using  the  measured  peak  current. 

u   =   ^  S^  W 

^       (1-R  )  qV    T  E  ^^^^ 

r  ^  appl 

where : 

V    1   =   applied  voltage 

4 .   Square-Pulse  Response 

Square-pulse  response  occurs  when  the  source  response 
is  longer  than  the  detector  carrier  relaxation  time.   It  is  a 
measure  of  the  relative  change  in  the  electron  and  hole  popula- 
tions.  The  transient  response  for  InP:Fe  is  controlled  by  the 
trap-assisted  recombination  of  electrons  and  holes  [Ref.  1]. 

For  square-pulses  (long  radiation  pulse),  the  uncharged 
acceptor  states  become  negatively  charged  due  to  electron 
trapping.   This  is  a  very  fast  process  which  leads  to  saturation 
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by  the  negatively  charged  acceptor  states.   At  the  same  time, 
hole  trapping  by  the  negatively  charged  acceptor  states  is 
occurring  but  at  a  slower  rate  due  to  a  substantially  lower 
hole  mobility.   At  a  certain  point,  the  hole  population  ex- 
ceeds the  electron  population  and  thus  makes  the  dominant 
contribution  to  the  current.   When  both  electron  and  hole 
trapping  rates  equal  the  hole-electron  generation  rate, 
steady-state  is  achieved.   [Ref.  1]   This  is  entirely  dependent 
on  the  intensity  of  the  incident  square-pulse.   The  peak 
current  response  is  measured  by  the  following  formula. 

(1-R  )  qv  p  T 

T   =   = - (23) 

^  S  W  ^       ' 

Parameters  are: 

p   =   optical  power 

T   =   recombination  lifetime 

5 .   Photoconductive  Gain 

Photoconductive  gain  is  defined  as  the  ratio  of  the 
recombination  lifetime  of  the  charged  carriers  to  the  transit 
time  across  the  active  area  between  electrodes.   The  bulk 
picosecond  photoconductor  decouples  the  dependency  of  the  carrier 
recombination  lifetime  with  the  transit  time.   The  response 
speed  is  a  function  of  the  carrier  relaxation  time  within  the 
active  region  and  is  thus  independent  of  electrode  spacing 
and  bias.   A  necessary  condition  for  photoconductive  gain  to 
be  greater  than  unity  is  provided  by  the  injection  of  electrons 
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into  the  device  by  the  metallic  contacts.   This  replenishes 
those  electrons  leaving  the  exit  electrode.   For  each  non- 
reflected  incident  photon,  more  than  one  electron  is  delivered 
to  the  external  circuit.   Photoconductive  gains  as  high  as  5 
have  been  demonstrated  in  InP:Fe  photocnoductors  with  AuGe 
and  AuSn  contacts  annealed  at  450  C  [Ref .  22] . 


41 


III.   EXPERIMENTAL  PROCEDURES 

A.   RADIATION  FACILITY 

Irradiations  were  performed  using  the  Electron  Linear 
Accelerator  (LINAC)  at  the  U.S.  Naval  Postgraduate  School. 
The  LINAC  is  a  travel-wave  accelerator  consisting  of  a  disk- 
loaded  circular  waveguide  thirty  feet  long  divided  into  three 
ten- foot  sections.   A  series  of  three  klystron  amplifiers  are 
used  to  feed  RF  power  into  each  ten- foot  section.   Electrons 
are  initially  produced  by  an  electron  gun  and  are  injected  into 
one  end  of  the  sections  where  those  in  proper  phase  relation 
will  be  accelerated.   A  magnetic  quadrupole  doublet  located 
downstream  from  the  deflection  magnets  serves  to  focus  the 
beam  on  the  target. 

The  LINAC  operates  with  an  electron  beam  pulse  length  of 
1.0  micro-second  at  a  repetition  rate  of  60  pulses  per  second. 

Electrons  are  focused  on  a  target  positioned  in  a  target  chamber 

—  fi 
with  a  vacuum  of  10    Torr  maintained  by  a  diffusion  pump.   In 

this  experiment,  electrons  leave  the  accelerator  at  an  energy 

of  100  MeV  and  an  average  electron  current  of  1  micro-ampere. 

Electron  fluence  is  measured  by  secondary  emission  monitors 

(SEM) .   Two  SEM's  were  used  during  the  experiment.   A  smaller 

SEM  was  placed  inside  the  vacuum  chamber  in  front  of  the  target. 

A  larger  SEM  was  placed  outside  of  the  vauum  chamber  in  the  beam 

line  after  the  target.   Prior  to  irradiation  of  the  target,  the 

SEM's  were  calibrated.   The  larger  SEM  was  previously  calibrated 
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against  a  Faraday  cup  and  was  determined  to  have  a  6%  efficiency 
[Ref .  23] .   The  smaller  SEM  was  calibrated  against  the  larger 
SEM.   As  the  electrons  impact  the  SEM's,  secondary  electrons 
charge  a  capacitor.   A  voltage  integrator  measures  a  voltage 
developed  across  the  capacitor  which  is  proportional  to  the 
charge.   The  ratio  of  voltages  between  SEM's  is  2.33  which 
equates  to  a  2.6%  efficiency  for  the  smaller  SEM.   The  electron 
fluence  is  determined  by  the  formula 


( . 026)  q A 

2 

where  $  is  the  electron  fluence  in  electrons/cm  ,  C  is  the 

capacitance  in  Farads,  V  is  the  integrated  voltage  in  Volts, 

q  is  the  charge  of  an  electron,  A  is  the  area  of  the  beam  in 

2 
cm  and  .026  the  SEM  efficiency  correction  factor.   A  more  ex- 
tensive description  of  the  LINAC  and  its  operating  parameters 
may  be  found  in  Reference  24. 


B.   SAMPLE  FABRICATION 

Devices  are  fabricated  from  two  different  types  of  materials 
Indium  Phosphide  with  Iron  dopant  (InP:Fe)  and  Gallium  Arsenide 
with  Chromium  dopant  (GaAs:Cr).   The  procedure  was  developed  and 
performed  by  the  Los  Alamos  National  Laboratory  (LANL)  Group 
E-11.   The  single  crystal  materials  are  obtained  commercially, 
the  InP:Fe  from  CrystaComm  Inc.  and  the  GaAsrCr  from  Cambridge 

Instruments  LTD.   The  InP:Fe  samples  are  claimed  by  the  manufac- 

7 
turer  to  have  a  resistivity  of  2  x 10   ohm-cm  and  an  electron 

2 
mobility  of  2200  cm  /V-sec. 
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The  crystal  is  cut  with  a  diamond  saw  to  a  specified  size 
(5  iranL  xi  mmH  x  .5  imnD)  .   The  surface  is  polished  and  bromine 
etched  to  remove  any  damage.   The  samples  are  then  rinsed  with 
methanol  and  de-ionized  water  to  remove  the  bromine  solution. 
Positive  resist  photolithographic  liftoff  is  applied  while 
spinning  (§   1000  rpm/5  sec  and  then  baked  for  30  minutes  at 
60°C.   The  samples  are  then  rinsed  with  methanol  while  spinning, 
baked  again  and  placed  in  a  masked  aligner  where  they  are  ex- 
posed to  ultra-violet  light  for  25  seconds.   Metallic  contacts 
are  applied  to  the  device  by  a  multi-step  vacuum  evaporation 
of  Au,  Ge,  Ni  compount  (80%  10%  10%).   Acetone  is  used  to  dis- 
solve the  .5  mm  masked  gap  exposing  the  material.   The  samples 
are  then  annealed  at  4  50 °C  for  2  min.   The  sample  is  supported 
by  a  5  cm  x  5  cm  alumina  substrate  with  50  ohm  microstrip  trans- 
mission lines  and  SMA  coaxial  connectors.   A  section  of  the 
transmission  line  is  removed  and  the  sample  is  mounted  with 
conductive  silver  epoxy  which  makes  electrical  contact.   During 
all  tests,  one  end  of  the  sample  connection  was  biased  and  the 
opposite  end  (current-signal  measurement)  was  connected  into 
the  sampling  oscilloscope  (25-ps  risetime) .   The  mounted  sample 
is  illustrated  in  Figure  8. 

C.   MEASUREMENT  TECHNIQUES  AND  INSTRUMENTATION 
1.   Pre- irradiation 

Before  irradiation,  all  samples  were  evaluated  to  deter- 
mine current  vs  voltage  characteristics  (dark  current) ,  impulse 
response  and  square-pulse  response  measurements.   A  total  of 
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Figure    8.      Mounted   photoconductive   detector    in    the 
double-ended    configuration 
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eight  InP:Fe  and  four  GaAs:Cr  samples  were  prepared.   Ingot 
numbers  are  listed  in  Table  3.   All  measurements  were  performed 
at  LANL. 

Dark  current  curves  were  determine  using  an  Ortec  Detec- 
tor Control  Unit  (Model  210).   The  InP:Fe  samples  were  biased 
up  to  500  volts  in  100  volt  increments.   The  GaAs:Cr  samples 
were  biased  up  to  125  volts  in  25  volt  increments.   Figures  9 
through  13  display  the  reference  plots. 

Impulse  response  measurements  were  evaluated  using  the 
equipment  listed  in  Appendix  A.   The  coaxial,  double-ended 
mounted  samples  (with  separate  bias  and  signal  cable  connec- 
tions) were  pulsed  with  the  Hamamatsu  Picosecond  Light  Pulser 
(9  0-ps  FWHM  optical  pulses)  at  a  frequency  of  10  KHz  and  a  wave- 
length of  820  nm.   During  the  procedure,  a  bias  of  200  volts 
was  maintained  on  the  InP:Fe  samples  and  100  volts  on  the 
GaAs:Cr  samples.   The  signal  was  evaluated  using  the  Tektronix 
TDR/Sampler,  Processor,  Acquisition  Unit  and  Oscilloscope. 
The  determination  of  amplitude,  full-width  at  half  maximum  and 
baseline  tail  was  made  and  photographs  taken.   Results  are 
listed  in  Table  3. 

Square-pulse  response  measurements  were  evaluated  using 
the  equipment  listed  in  Appendix  B.   The  samples  were  pulsed 
with  the  Tektronix  PG  501  Pulse  Generator  which  produces  an 
optical  source  with  fast  risetime  and  falltime  (150-ps),  con- 
stant power  (3  mW)  and  continuously-variable-duration  optical 
pulses  at  a  wavelength  of  830  nm  [Ref .  1] .   As  in  the  impulse 
response  measurements,  the  InP:Fe  samples  were  biased  at  200 
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volts  and  the  GaAs:Cr  samples  were  biased  at  100  volts.   The 
signal  was  processed  using  the  Tektronix  7104  Oscilloscope  (with 
listed  supplemental  units)  and  was  sent  through  the  Princeton 
Applied  Research  Signal  Averager  before  being  recorded  on  the 
Hewlett  Packard  7044A  X-Y  Recorder.   Both  photographs  of  the 
oscilloscope  signal  and  the  final  X-Y  graphs  were  used  to 
determine  the  rise  time  constants  of  the  square-pulse  response. 
The  results  are  listed  in  Table  4. 
2 .   Post- Irradiation 

After  the  above  characterization  was  made,  the  samples 
were  irradiated.   Two  coaxial,  double-ended  mounted  samples 
plus  a  phosphor  target  plate  were  mounted  on  an  aluminum  ladder 

and  placed  in  the  vauum  chamber  of  the  LINAC  for  irradiation. 

—  fi 
A  vacuum  of  approximately  10    Torr  was  maintained.   The  samples 

were  connected  to  coaxial  cables  to  allow  for  a  bias  voltage 

and  current  measurement  (dark  current)  during  the  irradiation. 

The  beam  was  focused  on  the  phosphor  target  plate  and  the  area 

was  adjusted.   A  video  camera  monitored  the  samples  to  allow 

for  proper  positioning  of  the  beam.   The  smaller  SEM  was  used 

along  with  a  1  micro-F  capacitor  to  measure  the  integration  of 

the  voltage. 

The  samples  were  irradiated  obtaining  fluence  in  the 

range  of  10    to  10    electrons/cm  .   Formula  (24)  was  used  to 

determine  the  exact  fluence.   Dark  current  measurements  were 

performed  just  prior  to  irradiation  and  immediately  after  the 

irradiation.   An  attempt  was  made  to  measure  dark  current  during 
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the  irradiation  by  placing  the  beam  in  standby,  however,  this 
was  fruitless  due  to  the  excitation  of  the  samples. 

After  the  samples  have  been  removed  from  the  vacuum 
chamber  and  monitored  for  residual  radiation,  they  were  sent 
to  LANL  for  post-irradiation  measurements.   These  measurements 
followed  the  exact  procedure  as  for  the  pre-irradiation 
measurements  and  the  two  were  compared  to  determine  the  effect 
of  the  electron  bombardment.   The  results  are  displayed  in 
Figures  9  through  25  and  listed  in  Tables  3  and  4. 
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Figure  9.   Dark  Current — InP:Fe   Ingot  number:   CCIPS-3148B 
Data  for  curve  #8  was  not  obtained  due 
to  device  breakdown  at  5  volts 
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Figure  10.   Dark  Current — InPrFe   Ingot  number:   CCIPS-2223 
Data  for  curve  #5  was  not  obtained  due  to 
device  breakdown  at  50  volts 
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Figure  11. 


Dark  Current — InP:Fe   Ingot  number:   CCIPS-3149T 
Data  for  curve  #6  was  not  obtained  due 
to  device  breakdown  at  39  volts 
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Figure  12.   Dark  Current--GaAs :Cr   Ingot  number:   Cambridge- 
A852/R 


52 


DARK  CURRENT  -  GAASiCR 


LEGEND 

g     ^i  A852/R-2  REFERENCE 

o     M2  A652/R-4  REFERENCE 
A     #3  A852/R-2  FLUENCE  4.34X10**14 
+"J4  A85VR-2"FLUENCE  L64X10**i5 
i' ■  |5  A8527R-2"FLUENCE  4.78X10**15 
o"  jir6"A8'52/R-4  FLUENCE '7.0rXiO**i5 


,f 


,A'- 


.-K 


-X  -  - 


200  300 

VOLTAGE  (V) 
100  MEV  ELECTRON  IRRADIATION 


400 


500 


Figure  13.   Dark  Current — GaAs : Cr   Ingot  number: 
Cambridge-A8  52/R 
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IMPULSE  RESPONSE  MEASUREMENTS 
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Figure  14.   Impulse  response  measurement  for  InP:Fe 
Drift  velocity 
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Figure    15 


Impulse  response  measurement  for  GaAs : Cr 
Drift  velocity 
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Figure    16. 


Impulse   response   measurement    for    InPrFe 
Full-width   at   half   maximum    (response    speed) 
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Figure  17.   Impulse  response  measurement  for  GaAs:Cr 

Full-width  at  half  maximum  (response  speed) 


59 


a)  before  irradiation  impulse 
response  FWHM  220-ps 


b)  after  irradiation  impulse 
response  FWHM  140-ps 


c)  before  irradiation  square- 
pulse  laser  response 


d)  after  irradiation  square- 
pulse  laser  response 


Figure  18. 


Comparison  of  InP:Fe  impulse  and  square-pulse 
laser  response  before  and  after  electron 
irradiation  for  a  representative  sample 
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a)  before  irradiation  impulse 
response  FWHM  260-ps 


b)  after  irradiation  impulse 
response  FWHM  140-ps 


c)  before  irradiation  square- 
pulse  laser  response 


d)  after  irradiation  square- 
pulse  laser  response 


Figure  19 . 


Comparison  of  GaAs:Cr  impulse  and  square-pulse 
laser  response  before  and  after  electron  irradi- 
ation for  a  representative  sample 
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Figure    20 


Square-pulse    response  measurement  for    InFrFe 
Peak   current   vs    Risetime    (10%-9  0%) 
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Figure    21. 


Square-pulse  response  measurement  for  GaAstCr 
Peak  Current  vs  Risetime  (10%-90%) 
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Figure  22.   Square-pulse  response  measurement  for  InP:Fe 
Peak  current  vs  Fluence 
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Figure  23.   Square-pulse  response  measurement  for  GaAs:Cr 
Peak  current  vs  Fluence 
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Figure  24 


Square-pulse  response  measureraent  for  InP:Fe 
Risetime  (10%-90%)  vs  Fluence 
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Figure    25.       Square-pulse    response   raeasurement    for   GaAsrCr 
Risetime    (10%-90%)    vs    Fluence 
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IV.   EXPERIMENTAL  RESULTS 

A.   DARK  CURRENT 

Dark  current  measurements  taken  before  and  after  the  100 
MeV  electron  irradiation  were  compared  to  determine  the  result- 
ing change.   Figures  9  through  13  display  the  results.   The 
mode  of  preparation  and  impurity  content  of  the  material  results 
in  variation  of  the  recombination  lifetime  [Ref.  16].   Since 
decay  times  are  inversely  proportional  to  the  Fe  dopant  levels 
in  InP:Fe  [Ref.  25] ,  samples  were  segregated  depending  on  the 
source  Ingot  and  whether  the  samples  were  prepared  from  the 
top  or  the  bottom  of  a  given  Ingot.   Higher  dopant  levels  in 
the  bottom  of  an  Ingot  are  due  to  gravitational  forces.   Compari- 
son of ref erence  dark  current  measurements  justify  this  segre- 
gation.  All  the  GaAs:Cr  samples  were  prepared  from  the  same 
portion  of  the  same  Ingot. 

Before  irradiation,  the  InP:Fe  samples  were  biased  up  to 
500  volts  in  100  volt  increments  and  results  indicate  a 
linearity  through  500  volts.   The  GaAs:Cr  samples  were  biased 
up  to  125  volts  in  25  volt  increments.   They  lose  linearity  at 
approximately  100  volts.   Post-irradiation  measurements  indicate 
a  noticeable  change  in  the  dark  current.   The  InP:Fe  sample 
results  indicate  an  increase  in  dark  current  for  the  applied 
voltages.   The  amount  of  increase  depends  on  the  fluence  with 
higher  fluence  yielding  higher  dark  currents.   The  GaAsrCr 
samples  show  a  decrease  in  dark  current  with  higher  fluence. 
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B.   IMPULSE  RESPONSE 

Impulse  response  measurements  for  pre  and  post  10  0  MeV 
electron  irradiation  are  listed  in  Table  3.   Included  in  these 
measurements  are:   1)  amplitude — from  which  I  (peak  current) 
can  be  determined  thus  allowing  for  the  determination  of 
electron  mobility  and  electron  drift  velocity,  2)  full-width 
at  half  maximum  (FWHM) — from  which  carrier  relaxation  time 
(response  speed)  is  resolved,  and  3)  baseline  tail--which  gives 
an  indication  of  whether  the  detector  returns  to  its  reference 
point  after  having  been  excited  by  the  photon  source,  and  if  so, 
within  what  time  period. 

Both  the  InP:Fe  and  GaAs:Cr  samples  are  affected  by  the 
irradiations  with  the  more  pronounced  effects  resulting  from^ 
the  higher  fluence.   Figures  14  and  15  are  graphs  of  the  effec- 
tive carrier  drift  velocity  derived  from  Formulas  21  and  22.   In 
these  equations,  electrons  are  assumed  to  be  the  carrier  type 
due  to  the  much  higher  drift  velocities.   The  measured  values 
of  the  peak  current  response  are  used  along  with  the  following 
experimental  parameters:   1)  reflectivity  (R  )  is  30  percent, 
2)  the  contact  spacing  (S)  is  equal  to  the  .5  mm  gap,  3)  the 
energy  necessary  to  create  an  electron-hole  pair  is  assumed 
to  be  5  eV,  4)  the  applied  voltage  (V    ^ )  is  200  Volts  for 
InP:Fe  and  100  Volts  for  GaAs:Cr,  5)  q  is  the  electronic  charge 
and  6)  the  optical  energy  deposited  into  the  detector  by  the 
incident  photon  source  (E)  is  approximately  580  pJ.   The  effec- 
tive drift  velocity  is  then  calculated  by 
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v^   =   yE^  (25) 


where  E^  =  4  kV/cm  for  InP:Fe  and  2  kV/cm  for  GaAs:Cr.   Both 
figures  show  a  decrease  in  the  carrier  drift  velocity  with 
increasing  fluence.   This  suggests  that  electron  mobility  and 
electron  drift  velocity  are  inversely  proportional  to  the 
incident  electron  fluence. 

The  response  speed  of  the  photoconductive  detector  is 
determined  by  the  full-width  at  half  maximum  (FWHM)  measurement 
of  the  impulse  response.   The  trend  of  decreasing  response  time 
(decreasing  FWHM)  with  increased  fluence  is  summarized  for  both 
material  types  in  Figures  16  and  17.   These  results  are  predic- 
table due  to  the  increased  number  of  trapping  and  recombination 
centers  created  by  the  displacement  damage.   Figures  18  and  19 
display  a  representative  sample  of  the  impulse  response  before 
and  after  irradiation  for  InPrFe  and  GaAs:Cr  respectively. 

For  the  GaAsrCr  samples,  all  baseline  tails  were  reduced. 
The  InP:Fe  sample  results  indicate  that  there  was  either  no 
improvement  or  an  increase  in  the  baseline  tails. 

C.   SQUARE-PULSE  RESPONSE 

Square-pulse  response  measurements  for  pre  and  post  100 
MeV  electron  irradiation  are  listed  in  Table  4.   Measurements 
include:   1)  amplitude — from  which  I  (peak  current)  is  determined 
and  2)  risetime  (10%-90%)  which  gives  an  indication  of  the 
time  involved  for  saturation  of  the  neutral  acceptor  states  by 
electrons  and  the  simultaneous  recombination  of  the  negatively 
charged  acceptor  states  by  holes. 
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Figures  18  and  19  display  a  representative  sample  of  the 
square-pulse  measurement  before  and  after  irradiation  for  InP:Fe 
and  GaAszCr  respectively.   Both  material  types  show  square-pulse 
nonlinearities  for  the  pre-irradiated  samples.   After  irradia- 
tion, the  GaAsiCr  samples  display  significant  improvement 
in  accurately  following  the  longer  events.   The  quality  of 
improvement  is  commensurate  with  higher  fluence.  This  response 
is  not  found  in  the  InP:Fe  samples,  but  rather,  an  increase  in 
nonlinearity  for  higher  fluence.   The  dependency  of  peak  current 
on  risetime  or  fluence  was  not  established  for  either  material 
type  as  displayed  in  Figures  20  through  23.   Risetime  vs  fluence 

is  plotted  in  Figures  24  and  25.   These  graphs  suggest  a  pro- 

14 
portional  relationship  for  fluence  above  2.81  x 10    (electrons/ 

2 

cm  )  for  the  InP:Fe  samples  and  an  inversely  proportional 

15  2 

relationship  for  fluence  below  4.78  x 10    (electrons/cm  )  for 

the  GaAs:Cr  samples. 

D.   COMPUTER  SIMULATION 

A  computer  simulation  of  a  100  MeV  electron  beam  incident 
on  a  rectangular  parallelepiped  (RPP)  semiconductor  target  was 
performed  using  the  CRAY  computer  located  at  the  Los  Alamos 
National  Laboratory.   The  ACCEPT  code,  one  of  the  ITS:   The 
Integrated  TIGER  Series  of  Coupled  Electron/Photon  Monte  Carlo 
Transport  Codes,  was  used.   The  ACCEPT  input  program  along  with 
the  Cross  Section  Generations  and  the  Monte  Carlo  executions 
for  both  the  InP:Fe  and  GaAsrCr  samples  is  enclosed. 
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Quoting  from  the  description  of  the  code,  "ITS  is  a  Monte 
Carlo  solution  of  linear  time-integrated  coupled  electron/ 
photon  radiation  problems.   The  ETRAN  model  was  used  as  the 
basis.   It  combines  microscopic  photon  transport  with  a  macro- 
scopic random  walk  for  electron  transport.   The  ACCEPT  code 
is  a  general  three-dimensional  transport  code  that  uses  a 
combinational  geometry  scheme"  [Ref.  26]. 

The  results  in  graphic  form  (normalized  to  one  particle) 
for  both  electrons  and  photons  are  displayed  in  Figures  26 
through  31.   Both  materials  give  similar  results.   It  is  obvious 
from  the  Angular  Distribution  graphs  (Figures  26  and  27)  that 
the  overwhelming  majority  of  both  electrons  and  photons  occur 
at  an  angle  between  0  and  10  steradians.   This  would  indicate 
that  there  is  very  little  deflection  of  the  incident  electrons 
or  the  resulting  photons.   The  Energy  Spectra  of  the  escaping 
electrons  (Figure  28)  indicate  that  most  are  leaving  the  target 
with  energies  between  90-100  MeV,  therefore,  only  a  small  per- 
centage of  the  total  energy  is  actually  deposited  into  the 

-2 
device.   The  actual  energy  deposition  per  electron  is  1.762  x lo 

11  -2 

MeV  (1.476  X 10    RADS-InP  Coul)  for  InP:Fe  and  1.993  x 10    MeV 

(1.500  X 10    RADS-GaAs/Coul)  for  GaAsiCr.   The  charge  deposition 

-3  -3 

is  -2.059  X 10    for  InP:Fe  and  -1.993  x 10    for  GaAsrCr.   Most 

of  the  photons  are  in  the  low  energy  range  (Figure  2  9) .   The  Flux 

Distributions  -(Figures  30  and  31)  substantiate  these  results 

with  the  highest  flux  for  electrons  in  the  high  energy  regions 

whereas  the  highest  flux  for  photons  is  in  the  lower  energy 

regions.   The  results  are  not  surprising  given  the  thickness  of 
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Figure  26.   Angular  Distribution  (Electrons) 
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Figure  27.   Angular  Distribution  (Photons) 
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ENERGY  SPECTRA  (ELECTRONS) 
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Figure  28.   Energy  Spectra  (Electrons) 
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ENERGY  SPECTRA  (PHOTONS) 
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Figure  29.   Energy  Spectra  (Photons) 
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ELECTRON  FLUX  DISTRIBUTION 
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Figure  31.   Photon  Flux  Distribution 
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the  sample  (thin  sample)  and  the  energy  of  the  incident  elec- 
trons thus  legitimizing  the  use  of  thin  sample  approximations. 
Included  in  the  Monte  Carlo  executions  are  listings  of  the 
physical  options  based  on  10  batches  of  100,000  histories 
each.   Bremsstrahlung  represents  the  highest  number  generated, 
more  than  all  others  combined. 

E.  DISPLACEMENTS 

In  determining  the  number  of  displacements  from  formulas 
5-10,  an  equal  probability  of  collision  for  the  secondary  atoms 
was  assumed.   Since  these  devices  are  compounds,  a  range  of 
the  displacement  cross  section  and  therefore  a  range  of  the 
number  of  displacements  is  accounted  for  through  the  Threshold 
Displacement  Energies  in  Table  1.   The  Kinchin  and  Pease  model 
was  used  in  determining  the  average  number  of  displacements 
per  primary  recoil  of  energy  T -, .   Table  5  is  normalized  to  one 
particle. 

F.  COMPARISON  WITH  NETURON  DAMAGED  DEVICES 

Results  obtained  from  electron  irradiation  of  InP:Fe  and 
GaAs:Cr  concur  with  those  obtained  by  Wagner  et  al.  [Ref.  1] 
for  neutron  irradiated  samples.   This  substantiative  agreement 
holds  for  dark  current,  impulse  and  square-pulse  response 
measurements.   Both  material  types  have  extremely  fast  impulse 
response  with  differences  showing  up  in  the  dark  current  and 
square-pulse  responses.   The  GaAsiCr  samples  exhibit  those 
desirable  qualities  at  higher  fluence  which  include  lower  dark 
current  and  better  accuracy  in  following  the  longer  event  of 
the  square-pulse. 
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TABLE  5 


DISPLACEMENTS 


Material 


In 


Ga 


As 


N  tlattice  atoms/ (an  )         29  09  70  79 

°  2.97x10^    2.97x10^    4.42  xio^    4.42x10^^ 


Minimum  Displacement 
Cross  Section  (cm^) 


3.94  xlO 


-22 


1.95  xlO 


-22 


Maximum  Displacement 
Cross  Section  (cm^) 


1.49  XlO 


-21 


1.07  XlO 


-21 


r^inimum  number  of 

Displacenents 

(atoms/cm-^) 


11.70 


43.66 


Maximum  number  of 

Displacements 

(atoms/cm-^) 


44.25 


47.29 


Minimum  number  of 
Displacements  in 
the  sanple  (atons) 


.029 


,109 


Maximum  number  of 
Displacements  in 
the  sample  (atoms) 


.111 


118 


Number  of  Replace- 
ments (atons/cm^) 


24.47 


25.03 


23.99 


23.61 
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V.   CONCLUSIONS 

A  new  class  of  compact,  reliable,  sensitive  and  extremely 
high  speed  radiation  detectors  are  currently  being  developed 
at  the  Los  Alamos  National  Laboratory  Electronics -Division. 
They  utilized  neutron  irradiation  to  introduce  trapping  and 
recombination  centers  into  the  semiconductor  to  speed  up  hole 
trapping  for  better  characterization  of  long-event  pulses.   In 
this  work,  we  have  demonstrated  that  electron  irradiation  at 
100  MeV  energies  is  capable  of  qualitatively  and  quantitatively 
producing  the  same  desirable  effects. 

Dark  current  in  the  GaAsrCr  devices  decreases  for  increasing 
fluence  whereas  dark  current  in  the  InP:Fe  devices  increases 
for  the  same  range  of  fluence.   Impulse  response  controlled  by 
carrier  relaxation  has  been  improved  by  electron  irradiation, 
with  response  speeds  <  100-ps  having  been  achieved  by  both 
material  types.   Electron  mobility,  drift  velocity  and  response 
speed  have  been  shown  to  decrease  with  increasing  fluence  for 
both  material  types.   In  contrast  to  InP:Fe,  the  electron  damaged 
GaAsrCr  devices  show  significant  improvement  in  accurately 
following  the  longer  event  square-pulses  by  exhibiting  linearity. 
The  GaAsrCr  devices  fulfill  the  necessary  requirements  for 
detector  utilization  and  are  expected  to  replace  InP:Fe  as  the 
bulk-semiconductor  material  for  many  applications. 

Electron  irradiation  from  the  NFS  LINAC  has  been  shown  to  be 
an  alternative  source  for  the  introduction  of  trapping  and 
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recombination  centers  resulting  from  irradiation  damage.  The 
LINAC  provides  the  advantage  of  a  well  focused  beam,  variable 
energies  and  controlled  fluence. 

ITS:   The  Integrated  TIGER  Series  of  Coupled  Electron/Photon 
Monte  Carlo  Transport  Codes  provides  a  theoretical  means  for 
characterizing  an  electron  beam  incident  on  detectors  of  various 
geometric  shapes  and  sizes. 
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APPENDIX  A 
IMPULSE  RESPONSE  MEASUREMENT  EQUIPMENT  LIST 

EQUIPMENT 

1.  Hewlett-Packard  6515A  DC  Power  Supply 

2.  Hewlett-Packard  3465A  Digital  Multimeter 

3.  Ortec  Detector  Control  Unit  Model  210 

4.  Hamamatsu  C1308  Picosec  Light  Pulser — 

Maximum  Peak  Power  10  Watts 
Wavelength  820  nm 
Repetition  Rate  10  KHz 

5.  Newport  Research  Corporation  Magnetic  Base  Model  100, 
430,  360,  420 

6.  Tektronix  7704A  Oscilloscope 

7.  Tektronix  P7001  Processor 

8.  Tektronix  A7704  Acquisition  Unit  with 

1)  S-6  Sampling  Head  and 

2)  S-53  Trigger  Recognizer 

9.  Tektronix  TDR/Sampler 

10.  Tektronix  Amplifier 

11.  Tektronix  Oscilloscope  Camera  C-27 
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APPENDIX  B 
SQUARE-PULSE  RESPONSE  MEASUREMENT  EQUIPMENT  LIST 


EQUIPMENT 

1.  Exact  Model  508A  Log  Sweep  Function  Generator 

Ramp  Mode  -  TRIG 
Ramp  Time  -  100  msec 
Mode       -  TRIG  PULSE 
Range      -  100  Hz 

2.  Hewlett-Packard  1900A  Pulse  Generator 

Output    5  0  Ohms 
Volts  into  50  Ohms 

Amplitude  2-5 

Offset     0-2 
Width  100-lK  nsec 
Time  Interval  10-100  micro-sec 
Rate  2.5K-2  5K 

3.  Tektronix  EG  503  Function  Generator 

Trigout  2.5V  into  600  Ohms 
Function    -  Square  Pulse 
Frequency   -  200  Hz 

4.  Tektronix  PG  501  Pulse  Generator 
Trigout  >  IV  into  50  Ohms 

5.  Hewlett-Packard  6515A  DC  Power  Supply 

6.  Hewlett-Packard  3435A  Digital  Multimeter 

7.  Keithly  177  Microvolt  DMM 

8.  Tektronix  7104  Oscilloscope 

a)  7S11  Sampling  Unit  with  Type  S-4  Sampling 

b)  7T11  Sampling  Sweep  Unit  ■ 

c)  7B92A  Dual  Time  Base  I 

d)  Oscilloscope  Camera  C-27 

I 

9.  Princeton  Applied  Research  (EG+G)  Model  4202  Signal  Averager 

10.   Hewlett-Packard  7044A  X-Y  Recorder 
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11.  Intronics  Regulated  Power  Supply  Model  EPM  200/15TD 

12.  Micro  Controle  Newport  Research  Corporation 

Model  B-1 
Model  AC-1 
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APPENDIX  C 
INPUT  PROGRAM  FOR  ELECTRON/PHOTON  TRANSPORT  ACCEPT  CODE 

A.   CROSS  SECTION 

1  MATE-RIAL  IN  .788  P  .212  DENSITY  4.78  SUBSTEP  40 

2  MATERIAL  GA  .482  AS  .518  DENSITY  5.32  SUBSTEP  40 

3  TITLE 

4  100.0  MEV  CROSS  SECTIONS  FOR  PHOTOCONDUCTORS 

5  ENERGY  100.0 


B.   MONTE  CARLO 

1  ECHO  1 

2  TITLE 

3  ...  100.0  MEV  KEIPPER  TEST  PROGRAM 

4  ************************  * GEOMETRY*  *********************** 

5  GEOMETRY  1  0 

6  RPP  0.00  0.50  0.00  0.10  0.00  0.05 

7  SPH  0.00  0.00  0.00  2,0 

8  SPH  0.00  0.00  0.00  4.0 

9  END 

10  Zl  +  1 

11  Z2  +  2  -  1 

12  Z3  +  3  -  2 

13  END 

14  0.0025  33.5  267.9 

15  *MAT  ECUT  PTCZ 

16  1 

17  0 

18  0 

2^9  *************************  SOURCE  ************************** 

20  ELECTRONS 

21  ENERGY  100.0 

22  POSITION  0.0  0.0  0.0 

23  RADIUS  0.5 

24  DIRECTION  0.0 

25  CUTOFFS  0.05  .001 

25  ******************* ****ouTPUT  OPTIONS******************** 

2  7  ELECTRON-ESCAPE 

28  NBINE  10 

29  PHOTON-ESCAPE 

30  NBINE  10 

31  ELECTRON-FLUX  1  1 

32  NBINE  10 

3  3  PHOTON-FLUX  1  1 

34  BNINE  10 

35  ********************** *oTHER  OPTIONS********************* 

36  HISTORIES  100000 
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ECHO  • 

TITLE 

...100.0  MEV  KEIPPER  RTES"  PR03SAW 


■GEOMETRY. 


ge6metry  1  0 

RPP  0.00  0.50  0.00  0.10  0.00  0.05 

SPH  0.00  0.00  0.00  2.0 

SPH  0.00  0.00  0.00  4.0 

END 

Z1  •♦I 

Z2  +2  -1 

Z3  +3  -2 

END 

0.0025  33.5  267.9 

•MAT  ECUT  PTCZ 

1 

0 

o 


ELECTRONS 
ENERGY  100.0 
POSITION  0.0  0.0  0.0 
RADIUS  0.5 
DIRECTION  0.0 
CUTOFFS  0.05  .001 

ELECTRON-ESCAPE 
NBINE  10 
PHOTON-ESCAPE 
NBINE  10 

ELECTRON-FLUX  1  1 
NBINE  10 
PHOTON- FLUX  1  1 
NBINE  10 

HISTORIES  1000000 


■SOURCE- 


-OUTPUT    OPTIONS' 


•OTHER    OPTIONS" 
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